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ABSTRACT 

This report is an analysis of a dual-plane amplitude-sensing mono- 
pulse system which relates to the essentially inertialess displacement of 
the sum and difference patterns in space with respect to the antenna coor- 
dinates. The development is predicated on the concept of scattering and 
transfer matrices. It is shown that the dual-plane electronic boresight 
scanning (EBS) problem considered here can be regarded as lying in an 
eight-dimensional vector space in which the terminal quantities are anal- 
ogous to the basis vectors which span the space. The transmission and 
reception properties of the circuit are separated by a direct-sum decom- 
position and only the reception case is considered in detail for the angular 
resolution problem treated herein. This is an interim report pertaining 
to the formulation of the problem and does not treat any particular phys- 
ical situation or present graphs of the output channel signals. A second 
report has been initiated for this purpose. 
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A CHAIN MATRIX ANALYSIS OF AN ELECTRONIC BORESIGHT 
SCANNING TECHNIQUE FOR A DUAL-PLANE 
AMPLITUDE-SENSING MONOPULSE ANTENNA 


INTRODUCTION 

The primary objective of a monopulse data-acquisition antenna is to provide 
a maximum gain in the data channel by pointing the peak of the sum pattern 
accurately in the direction of arrival of the telemetry signal. This is accom- 
plished by positioning the antenna so as to null the difference signal. When the 
RF system is properly aligned, the difference signal null coincides with the 
maximum of the sum signal and with the mechanical center line of the parabola. 
Proper alignment is important since the servo-control system senses only the 
difference-pattern null and assumes that the sum-pattern maximum is coincident. 

The purpose of the work reported here is to provide a means of separating 
the problem of directing the RF boresight accurately and rapidly from the problem 
of mechanically moving the antenna structure. The proposed technique requires 
the development of a method of scanning the sum-and-difference patterns (pre- 
ferably electronically) throughout a small cone centered on the mechanical center 
line of the dish. An equivalent effect may be achieved by operating on the signals 
from the RF feeds to simulate a motion of the boresight. This boresight control 
is referred to as Electronic Boresight Scanning (EBS). Given such a control, a 
conventional servo system, having reduced accuracy and response requirements, 
could be used to move the antenna structure to keep the line of sight to the space- 
craft within the EBS cone. The mechanical system would follow the spacecraft 
with moderate accuracy while the EBS system produced a correcting vernier 
motion (real or simulated) of the boresight relative to the mechanical center line. 
A brief discussion of the antenna servo-control problem will indicate the main 
basis for interest in the proposed system. 

A major restriction in the design of the servo-control system is produced by 
the mechanical resonances associated with the antenna structure, gearing, and 
drive. Typically, the lowest frequency resonances and, therefore, those of most 
concern, occur at two to three cycles per second in 85-foot antennas and produce 
resonant rises of 6 to 18 decibels (db). Both the damping and the frequency of 
the resonance may be expected to change with antenna orientation, wind loading, 
etc. Adaptive techniques to track and suppress resonances are being investigated. 
Lacking such advanced and complex systems, however, stability requirements 
dictate that the servo bandwidth be limited to frequencies substantially below the 
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resonant frequency. Figure 1 gives the open-loop frequency response (Bode 
plot) of a typical Type II system having a resonant peak due to the load. For a 
satisfactory transient response, it is necessary that the peak be not less than 6 
db below the 0 db gain level. The geometry of the plot determines the maximum 
obtainable cross-over frequency, F c , which is a close approximation to the 
closed-loop bandwidth. F c is normally less than one cycle per second in 85-foot 
antennas and may be 0.4 cycles per second or less. The shaded area indicates 
the frequencies to which the servo system must respond in following a satellite 
and correcting for disturbances. A low value of F limits the gain at these 
frequencies and, therefore, increases the error resulting from satellite track 
dynamics. To a lesser extent, errors caused by wind gusts and mechanical 
systems friction torque variations are also increased. 



Figure 1-Bode Diagram 

Additional design restrictions are imposed by stiction in the gears, bearings, 
and hydraulic -drive subsystem, particularly the servo values; drive-motor speed 
range, maximum torque, and torque ripple; and other drive and load non-linearities. 
These undesirable characteristics increase dead zone and degrade low-speed 
performance. 
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As a result of the restrictions listed, servo-control systems designed for 
the accuracy and transient response necessary in space data-acquisition antennas 
are complex and demand high levels of performance from critically adjusted 
components. This leads to systems which require frequent retuning and which 
inherently have less than the desired reliability. The proposed system is intended 
to relieve this situation by permitting a relaxation in the servo-control system 
specifications. Simultaneously, improved pointing accuracy and low-speed per- 
formance appear possible. 

One configuration of the EBS servo-control system is given in Figure 2. As 
has been stated, equivalent results, as far as the control system is concerned, 
could be obtained by actually scanning the sum and difference patterns about the 
mechanical center line of the dish or by operating on the signals from the feeds 
to give an apparent scanning effect. Figure 2 is based upon the latter concept. 

The system consists of two loops which are independent within limits except 
that the inner, EBS, loop has as its reference, or center position, the output of 
the servo loop. The EBS circuit has the property of producing a certain apparent 
angular displacement of the boresight per unit of E' applied. The EBS loop forms 
an accurate, fast-response, electronic-control system in which the boresight 
error, E, may be made to approach zero while the boresight scanning signal, E', 
remains as a measure of the misalignment of the center line of the dish and the 
line of sight. The boresight scanning signal is applied also to the servo-loop 
controller where additional compensation, required by the mechanical servo, is 
accomplished. The servo loop positions the dish to null E'. Since the two loops 
seek the same null conditions and since the bandwidth (and, hence, the speed of 
response) of the EBS loop is made greater than that of the servo loop, there is 
no destructive interaction between the loops. A brief transfer-function analysis 
of the system and results from an analog computer simulation are given in 
Appendix I. 



Figure 2-EBS-Servo Control System 
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Benefits which would be derived from a successful EBS development include 
improved low-speed performance due to the absence of mechanical stiction in 
the EBS loop. An essentially infinite speed range of the boresight would thus be 
possible. The greater bandwidth obtainable in the EBS loop would result in a 
reduction in the pointing error caused by target dynamics and load disturbances. 
In addition, step-response time would be decreased, improving acquisition capa- 
bility. Error in the EBS loop due to receiver noise is increased by the greater 
bandwidth, but since this is normally a minor contributor to total system error, 
Ref. [l], a large net gain should result. Perhaps the most significant benefits 
would be in system reliability. The decreased gain and bandwidth permitted in 
the servo loop allow components to be operated more conservatively and reduce 
the sensitivity of system performance to changes in component characteristics. 
There is less excitation of structural and hydraulic -drive subsystem resonant 
modes and, therefore, fewer failures due to vibration and fatigue. 


BACKGROUND FOR AMPLITUDE-SENSING MONOPULSE 

The problem under consideration is that of accomplishing electronic bore- 
sight scanning for a dual-plane amplitude-sensing monopulse system. Specifically, 
this is an attempt to determine the feasibility of displacing the sum and difference 
patterns in space at three frequencies (136 me, 400 me, 1705 me), simultaneously 
and independently, for the Goddard STADAN 85 -foot parabolic reflector instal- 
lations with 


F ( focal length) _ Q ^4 
d (diameter) 

At this time, structural modification of the existing installation is precluded and 
a technique which necessitates a minimum amount of retrofitting is sought. The 
existence of the conventional secondary patterns is assumed throughout. The 
problem considered here is that of operating upon, or processing, the signals 
derived from the secondary beams in such a manner as to bring about the dis- 
placement of the sum and difference patterns. In the absence of any known 
method for accomplishing dual -plane electronic boresight scanning, under the 
given constraints, the possibility of extension of a single plane technique de- 
scribed by Rhodes, Ref. [2], pp. 62-68 is investigated. The radio-frequency 
circuit associated with the single-plane technique is shown as Fig. 3. 

Quantitative results for electronic boresight scanning in the 85-foot reflector 
are obtainable in a straightforward manner from a mathematical model under 
certain idealizations. Values pertaining to absolute field patterns, error slope, 
and the relationship between physically realizable electronic phase shift and the 
pattern displacement in space at each of the three operating frequencies are 
derived from the mathematical analysis with suitable parametric variations. 
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Figure 3-Electronic Boresight Scanning Circuit For Single-Plane Amplitude- 
Sensing Monopulse Antenna 


Note: The figure above is extracted from the text "Introduction to Monopulse" by D. R. Rhodes 
McGraw-Hill (1959) 











Proceeding in a systematic manner, the overall problem is separated at the 
outset into two parts: (1) the monopulse four-beam secondary-pattern cluster 
or driving function and (2) the electronic boresight scanning circuit which suc- 
cessively transforms the four simultaneous signals of the driving function. 

The monopulse cluster of four secondary radiation patterns is developed as 
follows for each operating frequency. Referring once again to Reference [2 ], 
p. 35, the radiation pattern function for an aperture with cosinusoidal amplitude 
distribution and constant phase distribution for single-plane amplitude-sensing 
is given by 


77 2 


•cos u 


E =- 


rr* 

T 


- u 


u = 77 — sin 6 


where 


( 1 ) 


lim cos u 1 

u ~' 7T/2 _ 77 2 / 4 - u 2 

Here E is the field strength, d is the aperture diameter, A. is the wavelength, 
and 9 is the spatial angle measured from the normal to the aperture. It is noted 
that this radiation pattern is usually associated with a linear aperture, however, 
for purposes of this report the pattern function given by Eq. (1) is a satisfactory 
representative for the physically realizable patterns obtained in practice. A 
detailed calculation for the radiation pattern obtained with a nearly cosinusoidal 
current amplitude distribution on a circular aperture is given in Appendix H, 
together with a graph showing the minor deviation from the assumed relationship. 
The variation between the two functions is very small, at least as far as the -6 
db level, and differs principally in the region of the sidelobe structure. It will 
be shown that the effective tracking cone excludes the regions of significant 
deviation between these functions. 


If either the function 



cos u 



u 
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or E - A 2 (u) is assumed to represent the on-axis radiation pattern in three 
dimensions (i.e. if rotational symmetry is assumed to exist) then the polarization 
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problem has been ignored and an idealization has been introduced. The actual 
patterns obtainable in a parabolic reflector with crossed dipoles, say, have ap- 
proximately the same shape as the assumed form. In fact, an additional sim- 
plification can be introduced without adverse effect here although rigor is sacri- 
ficed. It will be assumed that the secondary pattern cluster can be represented 
by a beam function with the form of Eq. (1), but suitably displaced in space by 
the beam squint angle B % , ref. [2] p. 46. The feeds are offset from the focal 
point, but lie in the focal plane, so that 


sin 6 = K — - (2) 

s F 

where AX is the feed displacement, F is the focal length of the parabolic re- 
flector, and K is the beam deviation factor, ref. [3], p. 488. The value of K is 
approximately equal to 0.88 here since F/d = 0.424. Since the feed displace- 
ment is small, & s =r AX/Ffor the present problem. Amplitude-sensing by off- 
focus feeds is never quite pure due to the finite separation of phase centers. 

The driving function will be based on the assumption of coincident phase centers, 
however, in the interests of simplicity. 

GENERATION OF THE SECONDARY PATTERNS 
FOR THE DUAL-PLANE CASE 

The generation of the secondary pattern four-beam clusters utilizes the 
Eulerian angle concept of coordinate transformation; ref. [4], pp. 1-7-109, 
ref. [5 1 pp. 95-106. At each frequency, the assumed on-axis pattern is effec- 
tively rotated by the amount < 9 s , about a common center and into the proper 
quadrant, to produce one of the patterns of the beam cluster. A tilted coordi- 
nate frame is associated with each beam pattern in each cluster and the Euler 
transformation is utilized to relate field intensities on a titled space such as 
( r', O', cp') to the coordinate frame (r, 0,<p) attached to the parabolic reflector. 
That is, four titled spaces are related to the reflector coordinates to produce 
the complete four -beam cluster at a single frequency. The complete transfor- 
mation, in full generality, has the following form in Cartesian coordinates 
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where each element of the column matrices is converted to the corresponding 
spherical coordinate expression. Field intensities can be derived for the pat- 
terns of the cluster by selecting an arbitrary field point on a hypothetical unit 
sphere (1, 6 , cp), evaluating the unprimed column matrix, and performing the 
indicated product Ax. Then the primed column matrix is determined, and an 
ordered pair {O' cp') is found. In this manner the simultaneous cluster intensities 
are found for any trajectory in {0, cp). The matrix A is different for each pattern 
of the beam cluster at some specified frequency so that A - A. and 


x 1 = Aj x 
x 11 - A 2 x 

X nl = A 3 * 


( 4 ) 


for beams i = 1, 2, 3, and 4 respectively. A complete description of the Euler 
transformation, a sketch showing coordinate conventions and orientation of the 
antenna, and a table of the Eulerian angles a, fi, y for the four-beam clusters 
at 136 me, 400 me, and 1705 me are found in the Appendix. Determination of the 
primed coordinates ( 6 1 , cp 1 ), ( 6 11 , cp 11 ), (0 m , cp Iir )and ( 0 IV , cp IV ) allows evaluation 
of the expression 


cos u. 
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for = 1, 2, 3, 4 


Since u A = sin#., 6 - 6 1 , 0 11 , 9 111 , 0 IV . 

PROPERTIES OF THE WEIGHTING CIRCUIT 

The weighting circuit is fundamental to the EBS technique under considera- 
tion, therefore its properties are discussed in some detail. Consider a weight- 
ing circuit composed of the two hybrid junctions with the scattering matrix of 
the following form, 
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and an ordinary phase shifter, connected as shown in Figure 2. 



Figure 4-Weighting Circuit 


Let a, = I , a 2 e ' > e = A, represent the sum channel and error channel, respec- 
tively, in a single plane amplitude-sensing monopulse system. The time varia- 
tion ei is factored out and a^ > 0 , a 2 > 0. Also 9 = 0 or 6 = 7 t only, which 
is consistent with the fact that the input signals from the sum and difference 
network are in either a cophase or antiphase relationship as the source of radi- 
ation is on one side or the other of the unscanned boresight. Since each input 
signal gives rise to two output signals at any output port, each of and b 4 is 
the complex sum of four signals. It is possible to count phase shifts and allow 
for power splitting in each hybrid junction to arrive at the magnitude and phase 
of b 3 and b 4 relative to the input quantities a, and a 2 e'> e . A transfer matrix 
calculation T eg = 1 T. . ^T. B T. . can also be employed as will be explained in 
detail in the section entitled TRANSFORMATIONS OF THE CHAIN MATRIX IN 
V g (c) of this report. In either case, the following results are obtained. 


2 b , = a , - 


- a x e’ J ^i + ) a 2 e~ J + ja 2 e' JW e 
j a 1 + j a^'^ 1 + a^e'* 6 e’^ 1 


->e -Mi 


Using the identities 1 - cos 0 = 2 sin 2 0/ 2, 1 + cos0 = 2 cos 2 0/2, and 
sin 0=2 sin 0/2 cos 0/2 it can be shown that: 

0j 0 X 0 X 0 1/2 

I b = a 2 sin 2 — + 2a, a. sin — cos cos 9 + a 2 cos 2 

ji 2 12 2 2 2 2 

•Ai 0i 

= a x sin— + a 2 cos— for <9 = 0 , a x > 0 , a 2 > 0 


a, sm — - a„ cos — 
1 2 2 2 


for 6 - 7T, a x > 0, a 2 > 0. 
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a l 

0 

cos 2 — - 2 
2 

. ^ 0! 

a, a„ sm — cos — 
12 2 2 

CO s 

0 + 

a ^ 

sin 2 

= 

a l 

'/'i 

cos — ~ 

2 2 

. *1 

sin — 
2 
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Cb 

ii 

o 

a ! 

© 

A 

a 2 

> 0 
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a l 

cos— + a 2 

. *1 

sm T 

for 

Cb 

n 

3 

a l 

© 

A 

a 2 

> 0 


t an(ar g b 3 ) 


tan(arg b 3 ) 


tan(arg b 4 ) 


t an(arg b 4 ) 


+1 ( . *1 
cos — la. sin — 

2V 1 2 


*1 


+ a 2 cos — 


'Pi ( 

sin — (a sin — + a o cos ~~ 
2 \ ‘ 2 2 2 


cos — a, sin — 

2 V 1 2 


*1 


a 2 cos 2 


. 'Pi . 'Pi 

sin — la, sin — 

2 V 1 2 




a 2 cos 2 


'Pi f 'Pi . 0j 

cos y ( a i cos y a 2 sin y 

077 Ti 

sin — ( a, cos a 2 sin — 

2 1 2 2 2 


0! / 'Pi . 0! 

cos — a, cos — + a„ sin — 
2 \ 1 2 2 2 

. 0i "7 ^ r^v 

sm — a, cos — + a sin — 
2 \ 1 2 2 2 


for 0 = 0, a t > 0, a 2 > 0 


for <9 = 77, a } > 0, a 2 > 0 


for 6 = 0, a 1 > 0, a 2 > 0 


f o r 6 - 77, > 0, a 2 — ^ 


Cancellation of common factors in the numerator and denominator is avoided 
to preserve the quadrant of the phase angles (arg b 3 ) and(arg b 4 ). From the fore- 
going, it follows that the weighting circuit effectively introduces a phase advance 
of (7T — )/2 radians for the sum channel b 4 if |aj > |a 2 | . Let tan $ = tan(arg b 4 ) 
cot (0j/2) h tan (n/2 -\p x /2) where -v/2 < 0j <tt/2 is the range of 0j . An in- 
vestigation of t an(a rg b 3 ) shows that the difference channel information has been 
advanced by an equal amount 0 , to within an ambiguity of rr radians since the 
error channel inherently possesses a plus-minus ambiguity. The phasor diagrams 
of Figure 5, illustrate the phase and amplitude relationships for a few cases. 
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° 2 =1,0=0 a, =2, a 2 =1, 0 =180° Ql = 2, a 2 = 1, 6 = 180 

90° <A, =-90° (A, = + 90° 

CONSERVATION OF ENERGY CONDITION: a* + a^ =b^ + 

$ =(t7-^ 1 )/2 = PHASE SHIFT OF DUE TO WEIGHTING CIRCUIT 

Figure 5— Phase and Amplitude Relationships in a Weighting Circuit 
Pulse I (1-17), Pulse II (8-9) 




The first seven phasor diagrams depict the situation for a single-plane 
amplitude-sensing monopulse system (or one channel of a dual-plane system) 
when a plane wave-front illuminates the antenna from a point in space such that 
a single pulse of radiated energy is responsible for a sum intensity a x = 2 and 
a difference intensity a 2 =1. The parameter can take on all values over its 
range; only seven are shown. The modified difference channel Am now provides 
feedback information to minimize b 3 . Assuming that can be varied rapidly, 
so that the target position is sensibly unchanged in the time interval during which 
the error channel null is attained, the value of 4>i (including its sign) becomes 
the measure of target position when b 3 = 0. 

The eighth and ninth phasor diagrams illustrate the output relationships 
when the received plane wavefront arrives from a direction such that a x = 2, 
but a 2 = -1. These diagrams, together with the preceding seven, illustrate the 
nature of the v radian ambiguity associated with the phase angle of the error 
channel previously. Figure 5 shows that: 

1. Cophased input signals may result in cophased or antiphased output 
signals, 

2. Antiphased input signals may result in cophased or antiphased output 
signals. 

Equivalently, given b 3 = 0 for some 4 ) 1 (— tt/2 < <~p 1 < vr/ 2), there exist modified 
error signals (Am) of opposite sign on either side of the displaced boresight indi- 
cation for unmodified error signals ( A) of either sign. 

It is noted that an ideal or lossless weighting circuit necessarily yields a 
maximum modified sum intensity 2 = b . in terms of the total energy due to a 
and a 2 for a particular spatial coordinate 6 of the target when \p l is of such 
a value that b 3 = 0. In the event that the target moves to some other point & 2 on 
the hypothetical unit sphere, the density of the incoming plane wave remaining 
constant, the intercepts on the secondary patterns change, and the system energy 
takes on another value in general. As a consequence of this, </» x remaining con- 
stant, the maximum sum intensity will not always lie opposite the null indication. 
The phenomena described will be referred to as a separation of the difference 
pattern null and the sum pattern maximum is spatial coordinates (6, cp). The ef- 
fect is seen in the pattern plots of Figure 3, and the study of this effect for the 
dual-plane case is one of the objectives of the present analysis. In addition, the 
relative energy levels of all channels will be determined to provide a direct com- 
parison with the performance of an unscanned system. 


DUAL-PLANE ELECTRONIC BORESIGHT SCANNING CIRCUIT 

This report is concerned primarily with the analysis of one circuit and 
represents an attempt at extending the single-plane amplitude-sensing EBS cir- 
cuit described by Rhodes in Ref. [2] to a dual -plane EBS circuit. Figure 6. The 
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- 77/2 < \p 2 < 77/'2 


• T 1 3-T, # I 2 pT, y J 2 ,T 3 5 T, ^3 , 1 , 2 ,T 2 /3T3 0 T 4 ,T 4 / 4 

(a) MATHEMATICAL MODEL 



Figure 6 -Electronic Boresight Scanning Circuit For Dual-Plane 
Amplitude-Sensing Monopulse Antenna 


mathematical formalism described herein was developed for the purpose of in- 
vestigating alternative forms of EBS and provides a specific description of the 
function of every component in the circuit. It is possible to proceed first in a 
purely intuitive manner, however, by taking into account the basic properties of 
ring and square hybrids and, particularly, the properties of combinations of 
junctions such as the weighting circuit. For the case of two error signals, two 
weighting circuits appear to be a logical choice for displacing the bore sight null 
in each of two orthogonal planes in space. Splitting of the sum signal follows to 
obtain two weighting circuits which are roughly analogous to the single-plane 
circuit. The remainder of the dual-plane circuit is associated with phase nor- 
malization to compensate for (1) the phase changes introduced by the hybrid 
junction which splits the sum signal, and (2) the continuously variable phase 
changes introduced by the weighting circuits themselves, and finally (3) the 
phase change incurred at the combining hybrid. The latter should produce a co- 
herent sum channel similar to that of an unscanned system. 
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Junctions which exhibit the basic properties of ring hybrids and square hy- 
brids can then be arranged with suitable phase shifters and fixed delay lines as 
shown in Figure 6. The circuit of Figure 6(a) may be regarded as a mathematical 
model which ignores such issues as the minimum number of line crossovers, 
minimum line length, and physically realizable phase delays </> i > 0. It is noted 
that the unscanned or rest condition of this EBS circuit is given by \p lr - ^2 r = O.and 
-tt/2 < 4 > 1 £ tt/2, -tt/2 < <tt/2 is allowed at present in the mathematical 

model. A suitable physical version is presented as Figure 6(b), but is not utilized 
in the analysis. 


SCATTERING- AND TRANSFER-MATRIX NOTATION 

It is convenient, although not essential, to treat circuit problems such as that 
of Figure 6 in terms of the scattering-matrix analysis of ferrite devices. The 
notation of Lax and Button of the Lincoln Laboratories, ref. [6], pp. 506-539, will 
be used for most of this report with certain exceptions evidenced in the detailed 
development to follow. Additional background on scattering matrix notation can 
be found in ref. [7], pp. 401-479; ref. [8], pp. 516-530; ref. [9], pp. 350-375 and 
ref. [10], Chapter IV. Some discretion is required in making comparisons 
between the references cited above and others found in the literature as there 
appears to be no preferred or accepted set of conventions, which sometimes 
leads to confusion. Since the scattering matrix does not lend itself to multiplica- 
tion it is desirable to convert it to a transfer matrix that admits non-commutative 
multiplication with other transfer matrices. The conventions adopted in this 
report will be made definite by means of a detailed conversion from the scattering 
matrix to the transfer matrix for a junction having the form of a ring -type coaxial 
hybrid. (See Appendix). All other transfer matrixes used herein are derivable 
in exactly the same manner and are also tabulated in the Appendix together with 
their associated scattering matrices. 

The individual junctions such as ring hybrids, square hybrids, etc. can be 
described completely in terms of transformations in a 4 -dimensional linear vector 
space, ref. [8], in which the terminal quantities (i.e. input and output at each 
port) may be considered as basis vectors. A scattering matrix transformation 
relates the output at a given port to the inputs at all ports and its elements are 
complex reflection and transmission coefficients for a chain of junctions. A 
transfer matrix transformation relates an output or an input on the "right" of a 
chain of junctions to the inputs and outputs on the "left" of that chain. Elements 
of the transfer matrix are complex transmission coefficients only . A chain 
consisting of one junction is admissible for either the scattering or transfer 
matrices. The respective transformations are given by b = Sa" and g = Th. It 
is necessary to consider the transformation as being defined over the field (F) 
of complex values (C). Since the circuit of Figure 6 contains pairs of hybrids, 
the immediate problem of analyzing that circuit lies in an 8-dimensional vector 


14 



space v (F) = V 8 (C). It turns out that eight basis vectors can be made sufficient 

to span the vector space for this problem even though the dissipative or quadra- 
polar load termination of the conventional comparator circuit indicates a non- 
zero kernel, ref. [8], p. 142-146, or a variation in the number of ports at this 
interface. Ordinarily other means are required to succeed in problems of this 
type, and may involve rectangular matrices, ref. [9], p. 1. hi this instance, the 
problem is relatively simple and, since only one transformation corresponding 
to a physical termination occurs within the chain of fifteen transformations, the 
difficulty is avoided by regarding the transmission coefficient r of that load 
as a quantity arbitrarily close to zero. The inverse of T is then well-defined 
and square matrices are retained for the entire problem. It will be shown 
that the above restriction on r can finally be removed for both the transmitting 
and receiving cases. 


TRANSFORMATIONS OF THE CHAIN MATRIX IN V (C) 

o 

Referring once again to Figure 6, the associative property of matrix multi- 
plication is now used to separate the problem into several parts for convenience 
in displaying the product of 15 transformations of the input signals ( a x , a 2 , a 3 , 
a 4 ) provided by the antenna. It is noted that the operation of matrix multiplication 
is not, in general, commutative and requires that the order T ec) = T n T n _ x • • • 

= Tj T l4 • • • -Tj be preserved under the conventions previously adopted (Appendix) 
for the transfer matrices. From a consideration of the fundamental properties 
of the weighting circuit and the other components, the circuit arrangement given 
in Figure 6 should provide two orthogonal error channels similar to those 
described for the single-plane case in Rhodes. The sum input to each weighting 
package is reduced by a factor of 1//2 in voltage, however, and the combined 
phase correction 0 = (0 - 0 t )/2 should make possible a continuous coherent 
recombination of the modified sum signals emanating from the weighting packages. 
The values 0 X and 0 2 should be capable of restoring the null condition for both 
modified error channels simultaneously. That is, A^ = A az = 0 is sought and 
the modified sum channel signal should be maximized. The chain multiplication 
7 = T 15 T 14 • • -Tj is carried out for completeness and to verify the overall 

result. 

The conventional sum and difference monopulse network is readily described 
by a transformation # T 2 CT T X # T X , where an assumed form only is taken for the 
ring hybrid elements as indicated in the Appendix. Actual physical forms vary 
in practice, usually by some phase shift associated with the reference planes or 
by a crossover of the input or output cables. More involved variations can be 
identified, especially when waveguide is employed, but these will not be considered 
here since the specific description of a particular junction can always be accommo- 
dated. It is sufficient to utilize the form # T. of the Appendix for the ring hybrid. 
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Then 



where ^T, is an elementary matrix of the second kind, ref. [8 ], p. 193. By using 
the relation g = Th for the result above it is easily verified that 

_2b s = aj + a 2 + a 3 + a 4 , (2, Sum channel) 

-2b g = (aj + a 2 ) - (a 3 + a 4 ), (Aj, error channel) 
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-2b ? = (a t + a 3 ) - (a 2 + a 4 ), (Q, quadrupolar load) 

-2b g = (a t + a 4 ) - (a 2 + a 3 ), (A 2 , error channel) 

The equation g = Th also yields four similar output equations associated with 
the transmitting case, but these are not written explicitly since the present EBS 
problem pertains only to receiving phenomena. 

Next, the matrix representatives of the quadrupolar load termination (pT^), 
the cable crossover { c T 2 ) i the sum-channel power divider ( # T 3 ), and the phase 
normalization lines ) are annexed to the product by forming # T 3 a J 2 pT x ) 
which is given by: 



where ^Tj and pTj are elementary matrices of the first kind. Most of the matrix 
products above are of an extremely simple nature and lead directly to the result 


(/3 T 1 * T 3 a T 2 fV <t T 1 “ ^/3 T 1 a T 2 f T 1 # T 2 cr T l * 
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From this product, given in matrix form as 
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it is easily shown, via g = Th, that the original received signals have been trans- 
formed so that 


( a ! + a 2 ) - (a 3 + a 4 ) , 



(A,) 

~p~ ( a » + a 2 + a 3 + 

/2 

when Pj 

= 0, 

(2/72) 

— zr ( a i + a 2 + a ^j + a 4 ^ 
/2 

when Tj 

zz 0, 

(2/72) 

1 ( a j + a 4 ) - ( a 2 + a 3^ • 



(A 2 ). 


It is seen that the four channels are either in phase or at most anti-phase which 
is consistent with (1) the fundamental assumption that the signals (a t , a 2 , a 3 , a 4 ) 
derived from the antenna are ideally co-phased in a pure amplitude sensing 
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system and (2) the monopulse postulate that the sensing ratio for a positive angle 
of arrival is the inverse of the ratio for an equal negative angle, ref. [2 ], pp. 
20-26. (Multiplicative and additive sensing ratios, r m (u) = P(u)/P(-u) and 
r a (u) = P 0 (u) / P e (u) have been defined. P(u) and P(- u) are related to secondary 
beam intensities. P 0 (u) and P e (u) are odd and even functions derived from the 
former.) Subsequent development of the EBS network will lead to a violation of 
the postulate given above due to a lack of symmetry about the boresight. The 
first postulate, which states that monopulse angle information appears in the form 
of a ratio, will be preserved. 

In order to obtain two weighting circuits which are identical in all respects, 
except for the parameters < p , \p 2 related to boresight displacement in two 
orthogonal directions, the outputs thus far obtained will be permuted by means 
of matrix T. below. 

cr 3 



Matrix ( < T 3 ^Tj # T 3 1 2 pTj # T 2 ,T ) is the same as the preceding result 

with the indicated row permutation 1— — 3 , 2—4 and is therefore not written out 
here. The result obtained after operating with o J 3 will appear in the final result 
as one of two factors which completely describe the EBS circuit of Figure 6. 

The second factor is determined as follows. 

Forming the transfer matrix for both weighting packages, ,T 1 ,T 2 , 
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For simplicity in subsequent operations let 

, -i0i , J'P i 

Zj = 1 - e z 5 = 1 - e 


i "i'/'i , j ^l 

= 1 + e 1 z 6 = l+ e 


1 " j ^2 
z, = 1 - e 1 


1 j ^2 
z_ = 1 - e * 


z 4 = 1 + e 


■i'P 2 


z 8 = 1 + e 


^2 


Then the factor ( ^T 4 # T 4 a T 4 ^T 3 ,T 2 ^T 2 b Tj ) is formed using the following 
matrices: 


e - ^ 1 

1 


1 ' 1 
1 

1 

| 


l 1 

1 

1 

e-rt | 


l 1 

| 

e ' r | 


4 1 

H 

1 

T - 

’ cr 4 ” 
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1 1 

1 


! 1 o- -l 

1 i i i i 

1 1 

| 


1 l i i i 

l 1 -- 1-0 , 

1 1 


l 1 1 

1 1 0 

1 — 


— 1 _ 


where 0 = (0 2 - </i 1 )/2 
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1 
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Upon completion of the indicated operations, the two factors which describe the 
entire EBS circuit enter into a product, (^T 4 ,T 4 ct J 4 ^T 3 b T 2 2 B T X ) 

^ T 3 /3 T i • T 3 a T 2 r T i » T 2 ct T i -V’ and yield the final result T eq given below. 
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DIRECT SUMS AND REDUCTION OF RESULTS 


Prior to writing out the functional relationships for the elements of T eq , 
a rearrangement of the input and output terms of the column vectors (matrices) 
leads to a simplification. 




— 




b 5 


T 11 T » T ,S T 17 1 


a l 

b 6 


^31 "*33 ^35 ^37 1 


a 2 
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b 7 
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_ b 4_ 


Now g ' = T ' h ' and 
& eq- 



Tj ® T 2 , a direct sum of matrices, 
ref. [11] p. 25, ref. [12] pp. 165-169. 


The decomposition of T' into a direct sum and the accompanying decompo- 
sition of V g (C) into two invariant subspaces clearly distinguishes between the 
transmitting and receiving modes of operation. Since only the latter is of 
interest here the elements of T x only will be given. It is possible to set r\ 
equal to zero in subsequent expressions since the inverse of r x (or T" 1 ) does 
not appear in any of the relationships, Then the final output signals are given 
by: 


- 4b s = a i T 11 + a 2 T 13 + a 3 T 1 5 + a 4 T 1 7 ’ < A 1 > 

“ 4b 6 = a x T 31 + a 2 T 33 + a 3 T 35 + a 4 T 3? , (2) 
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where 


• 4b ? - a t T S1 + a 2 T 53 + a 3 T s$ + a 4 T S7 , 
. 4b g = T ?1 + a 2 T ?3 + a 3 T 7S + a 4 T ?7 , 
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T 7! =T„ + jZ 4 


T 73= T 75=7r- jZ <- 


It is of some interest to examine the final output signals for the case when 
the EBS circuit is in the rest condition corresponding to 0j = 0 2 = 0. Then 
0 = \p _ i/ j/2 = 0 also and, if a plane wave is incident on the aperture from 
the undisplaced bo resight, aj = a 2 = a 3 = a 4 . For a given beam cluster let this 
intensity equal a 0 . It is evident that 


b 5 = b 7 = b g = 0 (Aj, load, A 2 ), 

-4b g = 2 j (aj + a 2 + a 3 + a 4 ) =' 2 j (4a 0 ) (2) . 

From the latter equation, b g = -2j a 0 and the system energy Sj appears at 
the sum port (6) and equals |bJ 2 or 4aq. For this special case all of the 
energy appears at the sum port and is identical to that of a conventional mono- 
pulse installation with an ideal circuit. The value 8% = 4a^ will therefore be 
used as a reference in comparing the signal strength of the sum or telemetry 
channel for less restrictive examples (0 / 0, 0 2 ? 0, 0 i- 0 2 , source of 
radiation displaced from 6 - 0) since it represents the energy available when 
a mechanically scanned antenna is optimally directed toward the source of radia- 
tion in a tracking operation. 

Although the elements of Tj now completely describe the operation of the 
circuit and a computer could be utilized at this point to obtain field plots for 
arbitrary trajectories over the antenna, it is possible in this particular instance 
to achieve further simplification of the four expressions for the output quantities 

(F = sum, A 2 = error channel, 

A. = error channel, L~ load channel). 

It may not be desirable, or even possible, to reduce a given problem to the 

T eq . = T ® T 2 form in some instances, in which case the matrix product T n T n 1 

Tj can be solved by machine. 
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The simplification is effected for the receiving case by identifying factors 
related to the output signals from a conventional dual-plane sum-and-difference 
network. Using the direct sum decomposition for the receiving case, 



the following substitutions are made. 


2 

C 

72 


2/2 


(a x + a 


2 + a 3 


+ a 4 l 


A ic — - [(aj + a 2 ) - (a 3 + a 4 )] 
A 2c = - | [(a, + a 4 ) - (a, + a 3 )] . 


It is noted that the sign of A lc and A 2c depends only on the amplitudes of 
input signals (a.) associated with the secondary pattern cluster. 

Before proceeding with the simplification, it will be advantageous to review 
the character of the phase displacements due to a weighting circuit. Under 
-tt/2 - i//j < 7 r /2 and | a 3 | > |a 2 | , it can be seen that the common factors in 
the expressions for tan (arg b 4 ) always reduce to the quotient +1/ + 1 so that 
changes in the quadrant of (arg b 4 ) are due to the ratio cos ^^/sin ^ 1/2 • The 
numerator is an even function, therefore, the odd function in the 'denominator is 
solely responsible for changes in quadrant of $>. The values \p x = ±tt/2 deter- 
mine the bounds on the phase variation of \p as shown in Fig. 7. 
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4 > i = 0 




Figure 7-Sum-Channel Phase Diagram 


Figure 8-Difference-Channel Phase Diagram 


cos \p l/2 (aj cos \jj l/2 ± a 2 sin 0 1/2 ) 

tan (arg b 4 ) =— ; 7 7 : 7 r 

sin \p l/2 (aj cos </» 1/2 ± a 2 sin 0, /2 ) 

where aj > 0, a 2 £ 0 and should not be confused with the set a { of the input 

column vector. 

The expression for tan (arg b 3 ) is more involved due to the inherent 7 7- 
radian ambiguity. The ratios of the common factors can now be + 1 / + 1 or - 1 /- 1 . 
As in the preceding instance, the ratio cos ipi/2 /^i- n 0i/2 can change the 
quandrant of $. If the ratio of common factors is + 1 /+ 1 , then the range of $ 
is exactly as before. If the ratio is - 1 /- 1 , the region in the lower half of the 
phase diagram, Figure 8, contains the range of angles for the (arg b 3 ). 


COS 0 (a, sin 0 1/2 ± a 2 cos 0 ) 

tan (arg b,) — 

sin 1 p 1/2 (a J sin 0 1/2 ± a 2 cos 


where a t >0, a 2 ^ 0. 

The generalization to be made is that aj and a 2 will be interpreted as 
X c //2 and A 1C . A 2c respectively, and \p 1 / 2 may also become 0 2 / 2 , as 
required. 
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Then: 
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where 


$2 = tan -1 


COi *» /2 (71 


sin 0 2 /2 + A 2c cos0 2 /2 ^ 


sin 0 2 /2 sin i// 2 /2 + A 2 c cos 0 2 /2 ^ 


(5-t)" 

(I-t)- 


and A 2 c is P lus or minus > |^ c /v/2 I > |A 2c |. The £r (i = 1, 2, 3, 4) here 
represent the elements of the input column vector and are identical to the 
E. (i =1,2, 3, 4) in the development of the beam cluster in the Appendix. 
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where and $ 2 are the same as for the X channel. 
Since 


0 
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$( = ($ x ) or ($j + 7 t ) , 
$' = ($ 2 ) or C $2 + ' rT )< 


the following simplification is achieved. 
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It is noted that two signals are combined in a phase-coherent manner, above, 
to produce a sum channel (I) with unambiguous phase. That is, the 2 becomes 
a phase reference relative to A 1 , A 2 , L channels and may therefore be 
utilized for product detection or an analogous process. The A 2 and A x chan- 
nels possess the requisite 77 -radian ambiguity, as does the L channel, but the 
latter is of no significance in the present tracking operation. 


DISCUSSION OF RESULTS 

The results obtained in this report are as follows. An analytical description 
of a dual-plane amplitude-sensing monopulse EBS system has been obtained via 
scattering matrix notation and reduced to a simple form. Two satisfactory 
driving functions have been developed to represent the antenna under simplifying 
assumptions. The combination of driving functions and chain formulation make 
possible computer solutions for arbitrary trajectories of a satellite over an 
EBS-equipped monopulse installation. 
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The following can now be computed and presented in graphical form. 

(1) the exact relationship between the required amounts of electronic phase 
shift and the null displacement in spatial coordinates 

(2) the limits or extent of the null steering capability 

(3) all energy relationships among the sum and difference channels in terms 
of a comparable mechanically scanned system 

(4) the quantitative effects of null and sum-maximum separation 

(5) the effect on error-slope as a function of the scanning variables and 

'/V 


Preliminary antenna pattern cuts have been obtained via computer and these 
indicate that the schematic circuit is operating in a satisfactory manner. A 
second report, now in preparation, will exhibit numerous pattern cuts for the 
more significant trajectories over the antenna under various operating conditions. 
Fabrication of coaxial and stripline versions of the EBS circuit have been 
initiated. A signal simulator has been built to replace the antenna during the 
preliminary bench -test evaluations. A study of low-loss, high-speed electronic 
phase -shifters is proceeding concurrently with the work described above. Pre- 
sent plans include an antenna range evaluation of the complete EBS system 
utilizing a 15-foot parabolic reflector antenna. 


32 



APPENDIX I 


A brief transfer function analysis of the EBS-servo system is given in this 
appendix to support the heuristic argument presented in the introduction. In 
addition, an illustrative simulation of the system was made on an analog computer, 
and the results are reported here. Since the extent and the manner in which the 
servo specifications may be relaxed depend upon the outcome of the EBS develop- 
ment, a nominal relaxation, consisting of reductions in servo-loop gain and band- 
width to one half their normal values, was used. 


Transfer Function Analysis 


In Figure 9(a), the system of Figure 2 is redrawn to present the significant 
signals and transfer functions. The open-loop transfer function of the EBS loop 
is KGj while G 2 is that of the mechanical servo loop. K relates the pattern- 
scanning signal and the resulting scanning angle (S). For convenience K is 
assumed to be frequency independent and of unity magnitude. The effect of bore- 
sight sc anning is indicated by the summation in which the sum of the mechanical 
motion (M) and the EBS motion (S) produces the boresight (C). Figure 9(a) may 
be manipulated into either of the forms given in Figures 9(b) and 9(c). The 
relation between satellite line of sight (R) and the mechanical center line of the 
dish (M) is shown explicitly in Figure 9(b). If the gain of G, is much greater 
than 1 over the bandwidth of G 2 


and the closed-loop transfer function approaches that of the original system with- 
out the EBS loop, that is 


G 2 

ffc 

R 1 +G 2 


demonstrating that the EBS loop has negligible effect upon the mechanical servo 
loop. 


From Figure 9(c) the closed-loop transfer function relating the boresight 
angle (C) and the satellite line-of-sight angle (R) is 

c gi a + 

R " 1 + G t (1 + G 2 ) 
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SIGNA L 


(a) 



(b) 



Figure 9— EBS-Servo Control System Block Diagrams 
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by imposing the additional condition, |G 2 | «1 at frequencies near the mechanical 
resonances, 


R 1 +G X 


and system performance is essentially independent of the resonant effects. The 
conditions on Gj and G 2 are easily met since high-frequency correction com- 
ponents are carried by the EBS loop and the bandwidth of the mechanical servo 
loop is intentionally restricted. 


Simulation Results 


The Rosman, North Carolina, 85-foot antenna system was chosen for simula- 
tion as the servo loop in the illustrative example. Three system configurations 
were tested: 

1. The present Rosman system 

2. The Rosman system with reduced gain and bandwidth 

3. The Rosman system with reduced gain and bandwidth and with an EBS 
loop added. 

The closed-loop bandwidths were 1.6 cps for the normal Rosman system, 

0.82 cps for the relaxed Rosman system, and 5.0 cps for the EBS loop. The 
three configurations were tested for step function response, Figure 10; error 
in following a 0.01 cycle-per-second sine wave, Figure 11; and error due to a 
step-load disturbance, Figure 12. In each of these figures, the forcing function 
is recorded at the top, and the responses for the three system configurations 
are shown in order below. 

In Figure 10, the reduced gain and bandwidth are seen to produce an increase 
in settling time and a reduction in overshoot. Addition of the EBS loop restores 
the settling time and further reduces the overshoot. In Figure 11, the increased 
following error, resulting from the system relaxation, is shown while the bore- 
sight error in the system with EBS is almost imperceptible. The 0.01 cycle- 
per-second sine wave approximates the tracking requirements of a 100 -mile- 
high circular orbit. Figure 12 illustrates the reduced response to a load 
disturbance, such as a wind gust, when the EBS is added. Analysis of additional 
simulation data demonstrated reductions in hydraulic pressure surges and 
structural vibration. These effects may be inferred from the more sluggish 
response to the step input in Figure 10, as well. 
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APPENDIX H 

CIRCULAR APERTURE CALCULATIONS 


Following the development for the rectangular aperture of dimension a, b 
in Reference [ 3 ], p. 180 , the expression 

y*+a/2 j*+b/2 

g( 9 , 9) = F(^,T7)e^ k sinl9 (£ costp + T) sin9) d^d7J. 

•'-a/2 •'“b/2 


goes over to 




g( 6 , 9) = T F(r', (p')e jk ■*"«<'' «'<P' cosC <P + r ' .in<p'.in<p) r » dr ' d<p' 

*V = 0 *V=0 

for the circular aperture if a = radius, k = 2 n/\, and r'dtp'dr'= differential area 
on the aperture. Then 


,<p)= f f F(r',<p')e jkr ' sinecos( ^'P'> r' dr' d<p' 


g(^ 


•'(jp =0 ^ r =0 


if p = 


r'/a, u = 277 3 51 — , and F(r', 9') = F(r') only, 


\ 


/*2 * 

r 1 

g(u, 9) = a 2 

F(p) pdp 

*V=o 



J upcos(Cp-Cp') d , 



Figure 13— Circular Aperture Coordinates 
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For approximately cosinusoidal amplitude illumination with uniform phase 
distribution, 


F(p) = cos j = (1 ~ P 2 )- 

Therefore, since 



,jurcos(Cp-9') dcp » _ 2 7 T J Q (ur), 


g(u, <p) = 2 77 a 2 



(1 - p 2 )p [ J 0 (up)] dp. 


From Ref. [13], p. 163 


♦ W“ x > Jp<“) 


It follows, upon setting p = 1, a = u, and x = p that 


g(u. 9) 

277 a 2 


■ f 

J p=0 


Pj 0 (up)dp~ P 3 J 0 ( u P ) dp 


P= 0 


1 r l 2 f 1 2 

I pj 0 (up)dp+ p 3 J 2 (up)dp-— PJ 

Jp = 0 Jp=o ^=° 


Again, from Ref. [13], p. 162 


dx 


C xP J p ( ax )] = ax p J (ax) 


j(up) dp. 
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Upon setting p - 1, p-2, p = 3 successively in the above. 


i(u. jP) 

2na 2 


J ' (UP) 'I + 4 J 3 (UP) IJ ~ 2 ~2 J 2 (" P )\' 0 


^ + J 3 ( U )> J a (“) 


= -T J 2 ( u > 


U 


Since J N (u) = S^l An (u) , Ref. [14] , p. 128; 
N ! 


for N = 2, J 2 (u) = ^i_ A (u), and 


2 

g( u .?> = g(u) = A,(u). 

2 2 


By normalizing the above, the two radiation pattern functions which are to be 
compared are 


( 77 2 / 4 ) cos u 

f (u) = 

( n 2 /4 - u 3 ) 


and 


g( u ) - A 2 (u) , which are shown in Figure 14. 
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APPENDIX HI 

GENERATION OF SECONDARY PATTERNS 
FOR THE DUAL PLANE CASE 


x'= Ax 


A 

12 "13 


A n Aj 

Ajj k 22 Ajj 

^31 ^32 ^33 



Figure 15— Feeds and Beams 


x = r sin 6 cos <p x' = r sin 6 ‘ cos <p* 

y = r sin 6 sin <p y' = r sin 6 ' sin <p* 

z = r cos 6 z' ~ r cos 9 ' 

(cos y cos a-cos /3 sin a siny)(cos y sin a + cos /3cos a sin y)(sin y sin/3) 

A = | (-siny cos a-cos /3 sin a cos y)( - sin y sin a + eos/3 cos a cosy) (cosy sin /3) 

(sin/3 sina) (-sin^cosa) (cos /3) 


Euler angles 

a determined by feed positions. 

/3 equal to squint angle 6 s . 

y equal to zero due to assumed 
rotational symmetry of beams. 

Figure 16— Coordinate Conventions — 
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sin 9 * cos cp # = A n sin 0 cos cp + A 12 sin 6 sin cp + A J3 cos 6 

sin O' sin cp # - A 21 sin 6 cos cp + A 22 sin 9 sin <p + A 23 cos 6 ! 

cos^' — A 3 1 sin«9coscp + A 32 sin0sincp + A 33 cos 6 > 


E i = - 


. U = — Sin 0 '; E, => 6 1 , E, 0 11 , E, => 0 m , E 4 => 0lV 


cos u 



appendix IV 

SCATTERING AND TRANSFER MATRIX FORMS IN V 4 (c) 


i => row, j column 


Sjj Tj • Symbol 

tt 

1 input port 

output port 



(b) 










Symbo 1 


0 0 


0 0 


0 1 


1 0 


0 1 


1 0 


0 0 


0 0 


0 0 1 1 0 


0 0 I 0 1 


T. . = 

O' 


1 0 I 0 0 


0 110 0 


0 0 | Tj 0 


o oio r, 


r, o io o 


o r i o o 


r i o 


o r - 1 


o o 


o o 


o o 


o o 


r 2 o 


Figure 18-Scattering and Transfer Matrices 
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DERIVATION OF TRANSFER MATRIX FROM 
SCATTERING MATRIX IN V 4 (c) 


Given that 


S 


0 

0 

1 

1 


~ S n 

S !2 

S„ 

S 14 

0 

0 

1 

-1 


S 21 

s 22 

^23 

s 24 

1 

1 

1 0 
1 

0 


S 31 

^32 

S 33 

S 34 

1 

-1 

1 0 

0 


S 41 

S 42 

S 43 

S 44_ 


b = S a 


— 


— 

b l 


a l 

b 0 




= .s. . 

• 1J 

2 

b 3 


a 3 

_ b 4_ 


a 4 


h i - a i 


h 2 , bj 


h 3 , a 2 


h 4 , b 2 


1 


b 3 - gj 
a 3 ’ g 2 
b 4 ’ g 3 
a 4 ’ g 4 


Figure 19-Terminal Conventions 


/2 


b l " S 1X a i + S 12 a 2 + S 13 a 3 + S 14 a 4 


( 1 ) 


/2 

J" b 2 “ ®21 a i + ^22 a 2 + ^23 3 3 + ^24 3 4 


( 2 ) 


/2 


T~ b 3 ^31 3 1 + ^32 a 2 + ^33 3 3 + ®34 3 4 


( 3 ) 
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j b 4 - S 41 a i + S 42 a 2 + S 43 3 3 + S 44 3 4 ( 4 ) 

1 = Th 


g l 


h i 



T n T 12 1 
1 

T 

A 13 

T 

1 14 

g 2 


h 2 



1 

T 21 T 22 | 

T 

A 23 

T 

A 24 


= .T. 

• » j 



.T = 
• 1 j 

— - - 

- — 

— — 

g 3 


h 3 



T T 1 

A 31 A 32 

1 

T 

A 33 

T 

A 34 

g 4 


h 4 



^41 ^42 | 

T 43 

T 

A 44 






L ~ 


— 


g l 

= T u h . 

+ ^12 b 2 ' 

f T 13 h 3 + T 14 

h 4 



g 2 “ "^21 b l + ^22 b 2 + ^23 b 3 + ^24 b 4 (®) 

g 3 _ "^31 H X + ^32 b 2 + ^33 b 3 + ^34 b 4 (^) 

84 - T 4 i hj + T 42 h 2 + T 43 h 3 + T 44 h 4 (8) 

Solving equations (3) and (5) for g t : 

S 31 h l + S 32 h 3 + S 33 g 2 + S 34 g 4 = ( T ll h l + T 12 h 2 + T 13 h 3 + T 14 h 4> ~ j~ 

Solving equations (4) and (7) for g 3 : 

S 41 h l + S 42 h 3 + S 43 g 2 + S 44 g 4 = ( T 31 h l + T 32 h 2 + h 3 + T 34 h 4> 

Therefore T, 2 = T 14 = T 32 = T 34 = 0 

T n=S 3 i(j//2) = (j//2) = T 13 =T 31 

T 33 = S 4 2 (3/ / 2) = -(i/^2) 
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Solving equations (1) and (2) for g 2 and g 4 : 


/2 

j h 2 “ S 11 h l S 12 ^3 “ S 13 g 2 + S 14 g 4 


/ 2 

j— h 4 - S 21 hj - S 22 h 3 = S 23 g 2 + S 24 g 4 

It follows that 

f h 2 h 4 


g, =' 



/2 ^4 

?4= j V S 13 _S 23 



Therefore 


S,3 

_^23 

^14 

S 24 

' S ,4 

^24 

, S X3 

S 23 

T 

= T : 

1 23 

A 4 1 

T 

A 24 

CM 

If 


Tai ^1 + ^22 N + ^23 ^3 + ^24 N 


- T 41 + T 42 h 2 + ^43 ^*3 + ^44 ^4 


T 44 = (j//2) 


Then 


• ‘j 



"~1 

0 

1 

o’ 


0 

-1 

0 

-1 

j 





/2 

1 

0 

-1 

0 


0 

-1 

0 

1 
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